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Abstract 

A complete calculation of the pion-nucleon loops that contribute to the transition operator for 
NN —7- NNtt up-to-and-including next-to-next-to-leading order (N 2 LO) in chiral effective field 
theory near threshold is presented. The evaluation is based on the so-called momentum counting 
scheme, which takes into account the relatively large momentum of the initial nucleons inherent in 
pion-production reactions. We show that the significant cancellations between the loops found at 
next-to-leading order (NLO) in the earlier studies are also operative at N 2 LO. In particular, the 
l/rri]y corrections (with being the nucleon mass) to loop diagrams cancel at N 2 LO, as do the 
contributions of the pion loops involving the low-energy constants Cj, i = 1 . . . 4. In contrast to the 
NLO calculation however, the cancellation of loops at N 2 LO is incomplete, yielding a non-vanishing 
contribution to the transition amplitude. Together with the one-pion exchange tree-level operators, 
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the loop contributions provide the long-range part of the production operator. Finally, we discuss 
the phenomenological implications of these findings. In particular, we find that the amplitudes 
generated by the N 2 LO pion loops yield contributions comparable in size with the most important 
phenomenological heavy- meson exchange amplitudes. 
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I. INTRODUCTION 



The reaction NN — > NNir has been extensively studied both theoretically and exper- 
imentally over the past decades. However, the near-threshold regime is still not yet fully 
understood. After the first high-quality data for pp — > ppw [l| became available (further 
experimental data can be found in, e.g., the review article @], with the latest measurements 
in Refs. [1, 0]), it quickly became clear that the original models failed to reproduce the 
new data. For example, the model of Ref. 0] fell short by a factor of two for the reaction 
pp — > dir + and by an order of magnitude for pp — > pp-n . Various attempts were made to 
identify the phenomenological mechanisms responsible for this discrepancy. 

The first theoretical paper to explain quantitatively the cross section pp — > ppn° was 
Ref. 0]. The new contribution in p] originated from the short ranged, irreducible cur- 
rents constructed directly from the nucleon-nucleon potential. A phenomenological inter- 
pretation of this mechanism was provided in Ref. [zj , where the exchange of heavy mesons 
(mostly a and u) followed by a pion emission via a nucleon-antinucleon pair (the so-called z- 
mechanism) was calculated. The mechanism was also shown to provide the missing strength 
for pp — > dn + in Refs. jsl, 0|. An alternative mechanism is based on the pion-nucleon 
rescattering diagram where the off-shell pion-nucleon amplitude plays a crucial role. It is 



well-known that the isoscalar pion nucleon scattering length is very small — see Refs. [10 
for its most recent determination — as a result of a cancellation of individually sizable terms 
which have different energy dependences. It therefore appeared natural that in the off-shell 
kinematics relevant for the pion production reaction the amplitudes are significantly en- 
hanced. This mechanism was also shown to be capable of describing the experimental data 



in both pp — > pp7T° Hi, |l2| as well as pp — > dir + [13j reactions. At this point there was no 
way to decide which of the mechanisms described captures the correct physics. 

Since pion interactions are largely controlled by the chiral symmetry of the strong in- 
teraction, one might naturally expect that chiral perturbation theory (ChPT) provides the 
proper tool to resolve the above mentioned discrepancy. However, the use of the standard 
ChPT power counting, which is based on the assumption that all relevant momenta are 
effectively of the order of the pion mass, was not very successful. The first calculations in 
this framework were done at tree level up to N 2 LO for both pp — > pp-K° (13- 16l| as well as 
for pp — >• c?7r + 13|, Il7| . These studies revealed, in particular, that the discrepancy between 
theory and experiment increases for the neutral channel due to a destructive interference of 
the direct pion production and the isoscalar rescattering contributions at NLO in standard 



counting. In addition, some loop contributions at N 2 LO were found in Refs. 18|, |l9| to be 
larger than the NLO contribution, revealing a problem regarding the convergence of the 
standard ChPT power counting. 

It was soon realized that the large initial nucleon momentum at threshold p, p = \p\ = 
y/rriNTn^, which is significantly larger than the pion mass m n , requires the modification of 
the standard power counting. The corresponding expansion parameter in the new scheme is 

X = p/A x ~ 0.4, (1) 

with A x being the chiral symmetry breaking scale of the order of 1 GeV. Here and in what 
follows, this power counting will be referred to as the momentum counting scheme (MCS). 



This modification was proposed in Refs. [14J, |T7j while the proper way to treat this scale was 



first presented in Ref. [20] and implemented in Ref. (21|, see Ref. [2| for a review article. 



The MCS expansion is performed with two distinct parameters, namely the initial nucleon 
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momentum p and the pion mass m w , where m n /p ~ p/A x . The pion loop diagrams start 
to contribute at a given order in the expansion parameter, which can be identified based on 
the power counting, and, unlike the standard ChPT power counting, continue to contribute 
at all higher MCS orders. 

Due to the fact that the Delta-nucleon (A-N) mass splitting is numerically of the order 



of p, the Delta-isobar should be explicitly included as a dynamical degree of freedom 14]. 
This general argument was confirmed numerically in phenomenological calculations 22N24J] . 



see also Refs. [14], [17|, |21|, |25| where the effect of the A in NN — >■ NN-n was studied within 
chiral EFT. However, in this paper we focus on contributions from nucleons and pions only. 
The A degree of freedom will be included in a subsequent publication. 

In the MCS, pion p-waves are given by tree level diagrams up to N 2 LO and the cor- 



responding calculations of Refs. [20|, |26] showed a satisfactory agreement with the data. 
Meanwhile, for pion s-waves loop diagrams start to contribute individually already at NLO. 
However, they turned out to cancel completely both for the neutral j2l[ and charged j2?J pion 
production, a result which is reproduced in this paper. To obtain this result for charged pion 
production, it is crucial to consistently take into account a contribution related to nucleon 
recoil in the ttN vertex as explained in detail in Ref . 27| • As a by-product of the consistent 



treatment of nucleon recoil effects in Ref. [27] , the rescattering one-pion exchange amplitude 
at LO was found to be enhanced by a factor of 4/3 which was sufficient to overcome the 
apparent discrepancy with the data in the charged channel. The first attempts to study the 



sub leading loop contributions were taken in Refs. [28, 29 . 

In this paper we advance the analysis for NN — > NNn at threshold to N 2 LO. In partic- 
ular, we evaluate all loop contributions at N 2 LO that involve pion and nucleon degrees of 
freedom. A complete calculation of all operators at N 2 LO (tree level and loops) including the 
A degree of freedom, and the subsequent convolution with the pertinent NN interactions in 
the initial and final states will be reported elsewhere. We will show in this paper that also 
at N 2 LO significant cancellations occur and that only very few loop topologies contribute 
to the final amplitude. 

The paper is structured as follows. In Sec. [Ill we present our formalism and discuss the 
hierarchy of diagrams as follows from our power counting. The next two sections are devoted 
to a detailed discussion of the results for the loop topologies proportional to the axial- 
vector nucleon coupling constant qa to the third (Sec. IIII|) and the first power (Sec. HVl) . In 
particular, we reproduce the cancellation of all NLO terms found in Ref. [27| and demonstrate 
that a similar cancellation pattern also takes place among the loop contributions at N 2 LO. 
In the latter case, however, the cancellation is not complete. Sec. |V] contains a compact 
summary of the results of Sees. II I II and IIVI Here, we also give explicitly the finite loop 
contributions which survive the above mentioned cancellation. In Sec. IVII the regularization 
procedure for the loop integrals is outlined. In this section we also compare the finite 
pieces of our loops at N 2 LO to the size of the contact 4Nii operator estimated based on 
phenomenological calculations. Finally, in Sec. IVHI we summarize the results of the paper 
and discuss phenomenological implications of the observed cancellation of loop contributions. 
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II. FORMALISM AND POWER COUNTING 



A. Reaction amplitude and Lagrangian densities 

The most general form of the threshold amplitude for the pion production reaction 
Ni(p) + N 2 (—p) — > N + iV + 7r in the center-of-mass frame, can be written as: 

M th (NN -> NNir) = U x % [o x - a 2 ) ■ p + A 2 {B x X a 2 ) ■ pi (n + r 2 ) ■ c/>* 



+ (<?i + &2) ■ P 



Bii(n- t 2 ) + B 2 (t 1 xt 2 ) -0*, (2) 



where <7i j2 and ti j2 are the spin and isospin operators of nucleons 1 and 2. The final state 
pion's three-component isospin wave function is denoted by (f>, e.g. <fi = (0, 0, 1) for 7r°- 
production and <f> = (l,i,0) / \/2 for 7T + -production. 

However, as follows from the angular momentum conservation and the Pauli selection rule 
for the NN system, a final s-wave pion in NN — > NNir can be produced via two angular 
momentum transition channels only, namely 3 Po — ► l £>oS an d 3 Pi — > 3 SiS, where we use the 
spectroscopic notation 2S+1 Lj for the NN states while the lower case s corresponds to the 
I = pion partial wave in the overall cms. Therefore, the two spin-isospin structures in 
Eq. (T5]) are redundant, and the reaction amplitude, that acknowledges the Pauli principle, 



can be rewritten, without loss of generality, as [211, |30 



M th (NN -> iV jVvr) = X(?iX a 2 ) ■ p (n + r 2 ) ■(/>*+ B {a x + a 2 ) ■ p (r x X r 2 ) • 0* , (3) 

with A = A\ + A 2 and B = B\ + B 2 . To derive ([3]) we used the fact that the spin (isospin) 
matrix element for the operator O = z(c?i — a 2 ) ■ p (O — i (ri — r 2 ) ■</>*) is equal to that of 
O — (a± x a 2 ) ■ p {O = (ri x r 2 ) ■ </>*) for s-wave pion production. 

The amplitude A in Eq. contributes to 3 Po ^os, which is the relevant transition 
amplitude for neutral pion production in pp — > ppir . Conversely, the amplitude B in Eq. 
contributes to the charge pion production in pp — > dn + , driven by the 3 Pi — > 3 S\S transition 
operator. Furthermore, in some channels such as e.g. pp — > pnir + , both amplitudes A and 
B contribute in a certain linear combination. 

It is convenient to write down the threshold reaction amplitudes in the form where the 
relevant spin-angular structure of the initial and final nucleon pairs are shown explicitly 1 

M pp ^ pp7T o = 4:iA{S-p)X'\ 

M pp ^ d7T+ = -2V2iB (Sxp)-e. (4) 

Here e is the deuteron polarization vector, p is the unit vector of the initial relative momenta 
of two nucleons, and S = X 2 a y^Xil 'V%, and Z't = x\' a yX 2 >/V2 denote the normalized spin 
structures of the initial spin-triplet and final spin-singlet states, respectively. 

The main goal of this paper is to derive the contributions to A and B that originate from 
loop diagrams. The loop diagrams can be separated in two different kinds: the ones involving 
only pion and nucleon degrees of freedom and the ones involving A (1232) excitations in the 
intermediate states. In this paper we concentrate on the first kind only where we include all 



1 The connection of the amplitudes A and B to the observables is given in, e.g., Ref. [2 
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relevant contributions at orders NLO and N 2 LO in the MCS as detailed in the next section. 
An evaluation of all MCS N 2 LO operators containing explicitly the A will be presented in 
an upcoming publication. 

Our calculations are based on the effective chiral Lagrangian in which the lowest-order 
(LO) nN interaction terms read in cr-gauge 3ll . |32| (more details on the pion-nucleon La- 
grangian can be found, e.g., in Ref. |33|) 
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The next-higher order interaction terms have the form 
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In the equations above f n denotes the pion decay constant, is the axial- vector coupling 
of the nucleon, and N (it) corresponds to the nucleon (pion) field. The ellipses represent 
further terms which are not relevant for the present study. 

The Lagrangian density for the leading 4tt vertex also needed for the calculation reads in 
the a- gauge: 

The loop diagrams (at N 2 LO) lead to ultraviolet (UV) divergent integrals. These UV 
divergencies are removed by expressing the bare LECs accompanying the five-point contact 
vertex at the same order in terms of renormalized ones. As a consequence of the Pauli 
principle, only two independent linear combinations of Lagrangian contact terms contribute 
to the transition matrix elements A and B. We will denote the corresponding amplitudes 
by .Act and Bqt for future references. 



B. Diagrams and Power Counting 

In ChPT the expansion parameter is Q/A x where Q is identified either with a typical 
momentum of the process or m n . The key assumption for convergence of the theory is 
Q <C A x . As mentioned in the introduction, the reaction NN — y NNn at threshold involves 
momenta of "intermediate range" p ~ y/m^rriN larger than m n but still smaller than the 
A x ~ mj\r- In the MCS we are thus faced with a two-scale expansion. For near threshold 
s-wave pion production, the outgoing two-nucleon pair has a low relative three-momentum 
p' and appears therefore predominantly in S-wave. We therefore assign p' an order m n and 
introduce the expansion parameter 

P' rn„ p 
p p m N 



5 



LO: 



Direct 




Direct 



WT 



NLO: 



K 

/ \ 
/ \ 
/ \ 
-4 V- 

Type II 



\ / 
/ \ 
/ \ 
V- 

Type Ilia 



— ^ v- 

Type Illb 



T 

A-- 

/ \ 

-V- V- 

Typc IV 



Box a 



Box b 



FIG. 1: Complete set of diagrams up to NLO (in the A-less theory). Solid (dashed) lines denote 
nucleons (pions). Solid dots correspond to the leading vertices from and C^}, as given in 

(2) 

the main text, stands for the sub-leading vertices from whereas the blob indicates the 
possibility to have both leading and subleading vertices. The N N contact interaction is represented 
by the leading S-wave LECs Cs and Cy. The red square in the box diagrams indicates that the 
corresponding nucleon propagator cancels with parts of the ttN vertex and leads to the irreducible 
contribution, see text for further details. 



The diagrams containing only pion and nucleon degrees of freedom that contribute to 
the reaction NN — > NNtt up to NLO in our expansion, are shown in Fig. [TJ Details of 
the evaluations of each of the loop diagrams can be found in appendices \M and [B] The first 
two diagrams in the first line are sometimes called the "direct" one-nucleon diagrams in the 
literature, whereas the last (rightmost) diagram is called the rescattering diagram. We will 
discuss both next. 

At leading order one needs to deal with the "direct" pion emission from a single nu- 
cleon where the nucleon recoil irNN vertex of OJj ([6]) is necessary in order to produce an 
outgoing s-wave pion. In addition, at LO there is a rescattering operator with the Weinberg- 
Tomozawa (WT) tttcNN vertex which, however, contributes only to the charged pion channel 
due to its isovector nature. In order to clarify the counting in MCS, we will concentrate on 
the first of two "direct" diagrams in Fig. [TJ In this diagram each vertex attached to the 
pion propagator involves a momentum p. The pion propagator itself involves a momentum 
p, i.e. it is counted as p~ 2 , whereas the nucleon propagator only carries an energy oc m v . 
The m" 1 of the nucleon propagator cancels the factor m n of the s-wave pion production 
vertex, which counts as ~ m^p/rriN- Thus, counting the "momentum flow" in the vertices 
and propagators of the diagram, gives an order of magnitude estimate of the direct dia- 
grams (as well as the rescattering diagram) as p/m^. These diagrams are counted as LO in 
MCS. Traditionally the LO direct diagrams have been evaluated numerically by including 
the pion propagator in the distorted NN wave functions, i.e. only the one-nucleon-pion pro- 
duction vertex gives the transition operator. Numerically, in the traditional distorted wave 
Born approximation approach, the "direct" term appears to be significantly smaller than 
the estimate based on our naive MCS's dimensional analysis. This suppression comes from 
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FIG. 2: The irN — > ttN transition vertex: definition of kinematic variables as used in Eq. ([8j). 
two sources: first, there is the momentum mismatch between the initial and final distorted 



nucleon wave functions [2| - see also Ref. [34j for a more detailed discussion. Secondly, 



there are accidental cancellations from the final state interaction present in both channels, 
pp — > ppn° and pp — > dir + , that are not accounted for in the power counting. Specifi- 
cally, the NN phase shift in the 1 Sq partial wave relevant for pp — > ppn crosses zero at an 
energy close to the pion production threshold [35|]. All realistic iViV scattering potentials 
that reproduce this feature show in the half-off-shell amplitude at low energies a zero at 
off-shell momenta of a similar magnitude. The exact position of the zero varies between 
different models, such that the direct production amplitude turns out to be quite model 
dependent. The suppression mechanism of the direct term for the reaction pp — > dir + comes 
from a strong cancellation between the deuteron S-wave and D-wave components. Thus, it 
is not surprising that numerically the "direct" terms in both channels are about an order of 
magnitude smaller than the LO amplitude from the rescattering diagram, which turns out 
to be consistent with the dimensional analysis. Since this LO contribution is forbidden by 
selection rules for pp — > ppir while allowed for pp — > dii + , one understands directly why a 
theoretical understanding is a lot more difficult to achieve for the former reaction. 

At NLO, which corresponds to the order p 2 /m^, loop diagrams illustrated in Fig.rjjstart to 
contribute to the s-wave pion production amplitude. For the channel pp — > ppir the sum of 
NLO diagrams type II, III and IV in Fig. [T] is zero due to a cancellation between individual 



diagrams [21|. However, the same sum of diagrams II - IV gives a finite answer for the 
channel pp — > dn + [21]. As a result the net contribution of these diagrams depends linearly 
on the NN relative momentum which results in a large sensitivity to the short-distance NN 



wave functions [36]. This puzzle was solved in Ref. [27|, where it was demonstrated that for 
the deuteron channel there is an additional contribution at NLO, namely the box diagrams 
in Fig. rjj stemming from the time-dependence of the Weinberg-Tomozawa pion-nucleon 
vertex. To demonstrate this, we write the expression for the WT ttN — > ttN vertex in the 
notation of Fig. [2] as: 

xr j M I- (gg + 1) 

V-k-kNN — '0+ m 7r 

= 2m n + U-m^+E- + zO j - {e-^- + tO^j , (8) 

where we kept the leading WT vertex and its nucleon recoil correction, which are of the 
same order in the MCS, as explained below. For simplicity we omit the isospin dependence 
of the vertex. The first term in the last line is the WT-vertex for kinematics with the on- 
shell incoming and outgoing nucleons, the second term the inverse of the outgoing nucleon 
propagator while the third one is the inverse of the incoming nucleon propagator. Note that 
for on-shell incoming and outgoing nucleons, the expressions in brackets in vanish, and 
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FIG. 3: One-loop diagrams contributing to s-wave pion production at NLO and N LO. Notation 
is as in Fig. [TJ 

the ttN ttN transition vertex takes its on-shell value 2m 1T (even if the incoming pion is 
off-shell). This is in contrast to standard phenomenological treatments where Iq in the 
first line of (JSJ) is identified with m n /2, the energy transfer in the on-shell kinematics for 
NN — y NNtt, but the recoil terms in Eq. (jHJ) are not considered. However, p 2 /itln ~ rn n so 
that the recoil terms are to be kept in the vertices and in the nucleon propagator 2 . The MCS 
is explicitly designed to properly keep track of these recoil terms. A second consequence 
of Eq. (jSJ) is that only the first term leads to a reducible diagram when the rescattering 
diagram with the ttN —y ttN vertex is convoluted with NN wave functions. The second and 
third terms in Eq. (JSJ), however, lead to irreducible contributions, since one of the nucleon 
propagators is cancelled. This is illustrated by red squares on the nucleon propagators in 



the two box diagrams of Fig. |3j It was shown explicitly in Ref. [27| that those induced 
irreducible contributions cancel exactly the finite remainder of the NLO loops (II - IV) in 
the pp —y cLtt + channel. As a consequence, there are no contributions at NLO for both tc° 
and 7r + productions, see also our results in the two first rows of Tables fl] and [TTJ 

In this paper we extend the analysis of the previous studies and evaluate the contribu- 
tion from pion loops at N 2 LO. Once the complete calculation at N 2 LO is performed, the 
calculated theoretical uncertainty based on our power counting is going to be reduced to 
~ (m^/m/vr) 3 / 2 < 10% for the amplitudes. At N 2 LO, one gets contributions from two sets 
of loop diagrams which differ in the power of g&- The diagrams proportional to g\ are the 
subleading contributions to the NLO diagrams of Fig. [TJwe already discussed. In addition, 
there is a set of pion loop diagrams proportional to g&, see Fig. |3j A naive MCS estimates 
indicate that the diagrams proportional to g^ could play a role already at NLO. However, a 
more careful analysis reveals that the contributions of each of these g^ diagrams at NLO is 
zero, see Appendix IB II for a detailed discussion. In subsequent sections it will be shown that 
partial cancellations take place among g^ and g\ diagrams at N 2 LO. Unlike the cancellation 
among the g\ diagrams at NLO, the cancellations at N 2 LO are not complete so that there 
is a non-zero transition amplitude from the gA and g\ diagrams. 



How to deal with the nucleon propagator in the MCS was shown in Ref. 28 1. 
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FIG. 4: Exemplary irreducible diagrams that contribute at higher than N 2 LO. Notation is as in 

Fig.m 

We already discussed pion loop diagrams with the nN — > tcN vertex stemming from the 
leading Weinberg-Tomozawa term, and its recoil correction, C^. In addition, there 
are two kinds of loop diagrams in Fig. [3] which involve the vertices from C\Jj\ those where 
the Q terms appear at the vertex, where the outgoing on-shell pion is emitted, and those 
where they provide an intermediate interaction. The former kind appears to be suppressed 
for s-wave pion production due to the pion kinematics near threshold. The Lagrangian term 
containing C4 can not contribute at an outgoing s-wave pion vertex since it is proportional 
to the gradient of the pion fields, cf. Eq. (Q. The contribution of the C3 Lagrangian term 
via this type of vertex is only non-zero if the term proportional to the time derivative of the 
pion field is considered in the Lagrangian, cf. Eq. ([6]). This C3 term, however, results in a 
loop amplitude which is suppressed by m n /m N compared to the leading loop at NLO, and 
thus it is of higher order (N 3 LO). In addition, the contributions proportional to LECs C\ and 
c 2 are quadratic with and therefore strongly suppressed. 3 However, the contributions of 
the vertices proportional to C2, C3 and C4 are potentially important at N 2 LO once embedded 
in the off-shell intermediate irN vertices (on the lower nucleon line of the gu-type diagrams 
in Fig. [3]). In what follows we will discuss the individual contributions of loops in detail. 

As a final remark, we give in Fig. H] some examples of additional loop topologies which 
start to contribute at a higher order than what is considered in the present study. The 
common feature of these diagrams is the presence of only one pion propagator inside the 
loops. As a consequence, by using appropriate integration variables, one can eliminate the 
large initial three-momentum p from the loop integrals, meaning the loop momentum will 
scale with m n . This explains why these loop diagrams in Fig. H]only start to contribute at 
order N 3 LO or higher. 

III. CALCULATION OF DIAGRAMS PROPORTIONAL TO g\ 

The diagrams of the gj^-group in Fig. [3] have a common structure illustrated in Fig. 
The loop diagram in Fig. [5] is integrated over the momentum I = (Iq,1). We also use the 
short-hand notation 

I — I + ki — q , 

with ki = pi — p[. The pion isospin indices a, b, and c are defined as shown in Fig. [51 The 
circle containing the vertex operator A abc produces an outgoing pion of isospin index a off 



3 Naively, the term proportional C2 scales as m w (Zo+recoils). However, a similar mechanism as the one 
explained below Eq. (8) forces the vertex to become proportional to roj. 
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FIG. 5: The general structure of g\ diagrams. 

nucleoli 1. This operator is different for each diagram and its explicit form is derived in 
Appendix |AJ where also the detailed structure of each ^-diagram is given. 

The invariant amplitude for each relevant diagram proportional to g\ can be written as 

where B 2 (l, /) is the common operator structure associated with nucleon 2 in Fig. |5j The op- 
erator B2(l, I) involves two pion propagators, two 7rNN- vertices and the nucleon propagator. 
The explicit form of B 2 (l, I) can be read off from the diagram in Fig. 

B 2 {l,l) 




y ^(s 2 .l- S '<- i2 Pj-»-> ), (10) 

h V 2m N J 

where v M = (1, 0) is the nucleon four-velocity and 5 M = (0, a/2) is its spin- vector. Note that 
£>2(M) contains no isospin indices as all isospin operators are included in Eq. (J9j). Since 
the structure of -E?2(M) is the same for all considered g\ diagrams we will consider, we 
concentrate our discussion on the structure of operator A a ^f in Eq. ([9]), see Appendix |A] 

for details. Note that the amplitude, Eq. ([§]) is not yet properly symmetrized with respect 
to the two nucleons. Below we will first discuss, how the partial cancellation amongst the 
various pion loop diagrams emerges on the basis of the decomposition illustrated in Fig. [51 
In Sec. |V]the non- vanishing remainder will be given in a symmetrized form. 



A. Pion s-wave contributions oc g\ 

In Appendix [A] we derive the expressions for each of the six g\ diagrams which contribute 
to near-threshold s-wave pion production from two nucleons. The results of these calculations 
are summarized in Table [I] where, for convenience, we have introduced the following short- 
hand notation for the isospin structures: 

r + = (r 1 + r 2 ) a , r_ = (r 1 -r 2 ) a , r x = z(ti x r 2 ) a . (11) 
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The left column in Table |T] shows the spin structures that emerge in these diagrams, the 
next six columns represent the contributions from the individual diagrams to the given spin 
structure, whereas the last two columns summarize the net effect of all diagrams and the 
MCS order, respectively. When we add the resulting expressions for the six diagrams we 
confirm the finding of Ref. [23] that the sum of the NLO contributions from all diagrams 
vanishes, see the first two rows of operators in Table [H Moreover, since the sum of the 
operators in the first two rows of Tabled] is an exact zero, the corresponding spin- momentum 
structures S% ■ I and S% ■ I will not contribute also at N 2 LO and all higher orders. In addition, 
all nucleon recoil corrections oc l/(2mjv) to the individual diagrams at N 2 LO also cancel 
in the sum. The reason for that cancellation is completely analogous to the cancellation 
that happens at NLO, see discussion below Eq. (jSJ). In fact, only those parts of the g\ 
diagrams that cannot be reduced to the topology of the diagram II in Fig. [3j give a non-zero 
contribution to the transition amplitude. Thus, only very few N 2 LO contributions to the 
pion production amplitude remain, as seen in Table HI The non- vanishing terms appear from 
the two cross-box diagrams and diagram IV. 

Since the sum of the A a ^f operators from the different diagrams starts to contribute at 

N 2 LO, we keep only the leading part of the operator B 2 (l,l). Adding up the contributions 
from all six g\ diagrams we arrive at the following result: 



iM f LO = i- 



g% f dH S 2 -l S 2 -l 1 



« 4/ ff 5 J (2vr) 4 P - ml + tO p - m 2 + zO -v ■ I + iO 

x{(-2. + + r,)-^- Q ( Sl ■ I) + (-2r + - rj-^fr ■ I) 

-8rASi-h) J l+ Jl- q \, (12) 



k\ — m 2 + iO 

where for the nucleon propagator in Eq. f JTO]) we dropped p 2 o and all recoil terms of order 



Oim^ compared to the lower-order Iq = v ■ I ~ \l\ ~ p term. Rearranging the isospin 
structure we arrive at three independent integrals to be evaluated for s-wave pion production: 



*Mf LO = ■ q )r + S^j dH lM + lh 



(27r) 4 (/2 _ m 2 + iQ )([2 - m 2 + iO)(- V ■ I + iO) 2 



(27t) 4 (/2 _ m 2 + iQ )(/2 _ m 2 + i0 )(- V ■ I + iQ) 



2 



x _ S 2 S 2 (Si ■ ki) f dH IMl + Qx 



x kf-ml + iOj (2tt) 4 (/2 _ m 2 + , 0) (/ ~2 _ m 2 + , ) (_„./ + t 0)j 

Employing dimensional regularization and an integration method outlined in Appendix O 
Eq. (Til?]) can be brought into the more transparent form: 

*Mf LO = 9 -^^^r +l e^k 1 ,S lu v a S 2 , [-J{k\ 



+r x (Si • h) 



24 v u 9(4tt) 



(14) 



where we have only kept the lowest order parts of the integrals which give contributions to 
the amplitude at N 2 LO. The pion loop diagrams generate ultraviolet divergent terms, which 
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TABLE I: Interference pattern of NLO and N 2 LO s-wave contributions from the individual g\ 
diagrams. The Table shows the contributions to the vertex t^t\ A ^ defined in Eq.Q and Fig. 

9a 

These contributions are given separately for the different spin-momentum structures of the vertex 
A a s c , shown in the leftmost column. The notation for the isospin structures is defined in Eq. tllll) . 

Si ■ I 
Si •/ 

Sl ■ (Pl+Pl') 2^ 
* ■ (Pl+Pl') ^ 

Si ■ f- 
Si -fc 



-u-i + iO 

1 1>2 _. m 2 _J_ A 



Type II 




Type Ilia 


Type 1Kb 


Type IV 


Box 


a 


Box b 


Sum 


Order 


-4r + — 4t_ + 


2r x 







-2r + - r x 


6t+ + 6t_ 


-2t_ - 


- T X 








NLO, N 2 LO 


4t + + 4r_ + 


2r x 


2r+ - 


T X 





— 6t_|- — 6t_ 







2r_ - r x 





NLO, N 2 LO 


-2t + + 2r_ 




2t+ - 


T X 













-2r_ + r x 





N 2 LO 


2r+ - 2t- 









-2r + - t x 





2r_ - 


rT x 








N 2 LO 


4t + + 4t_ — 


2r x 







2r+ + r x 


— 6t+ — 6t_ 


2t_ - 


hTx 








N 2 LO 


-4t + — 4t_ — 


2r x 


-2r+ + 


Tx 





6-rf + 6t_ 







-2r_ + r x 





N 2 LO 












-2r + - r x 













-2t+ - T X 


N 2 LO 







-2r+ + 


Tx 
















-2t+ + T X 


N 2 LO 















-8r x 










-8t x 


N 2 LO 



are contained in the following integral: 

jm = _ % f I (15) 

K l> J {2TT)* I 2 - ml + tO {I + h) 2 - ml + iO 1 ' 

The divergences are to be absorbed by the LECs accompanying the five-point (4iV7r) vertices 
as we will discuss in Sec. lYU 



IV. CALCULATION OF DIAGRAMS PROPORTIONAL TO g A 

We evaluate the qa diagrams following a similar strategy we used when we evaluated the 
g\ diagrams. The invariant amplitude for each diagram proportional to g A can be written 

as 

dH 



^ M 9a = D 2 (l, l) e bcd 4 < c , (16) 



where Dzil, I) is a common operator structure which is associated with nucleon 2 in Fig. [6j 

a « 



abc\ „/ 



Pi — p'i 

b \ _ / c / 

P2 p' 2 



FIG. 6: The general structure of qa diagrams. 

This structure involves the WT vertex at the second nucleon and the two pion propaga- 
tors: 

m !) = P-K + i P-rL + iO V ~^T- (17) 
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Note that we have only written the leading WT- vertex contribution in D 2 (l, I), Eq. (fl7|) . 
since the sum of the A a ^ operators starts to contribute at N 2 LO only, as can be seen in 
Table [TTJ This is in full analogy to the sum of the A a \ c operators, which also only start to 

9 A 

contribute at N 2 LO, where, as discussed just before Eq. (Tl2|) . the recoil (l/raw) corrections 
in B2(l,l), Eq. f JT0|) . only contribute at higher order. In other words, the corrections to 
D 2 (l, I), that is the recoil correction to the leading WT interaction term and the correction 
stemming from the C4-vertex, contribute at a higher order than what is considered in this 
work. Notice further that the c 2 and C3-vertices in Eq. ([6]) are isoscalars and, therefore, do 
not contribute to the function D 2 (l, l). The contributions of these LECs will be discussed 
in the next section. 



A. Pion s-wave contributions cx gA 



The operator expressions for each individual diagram of the g^-type contributing to s- 
wave pion production can be found in Appendix [B] In a complete analogy with the g\ 
diagrams, we summarize in Table [TT] the contributions of the individual diagrams and their 
net effect for different spin structures. In distinction to the g^-graphs the diagrams of this 
topology do not appear, contrary to naive MCS expectations, at NLO, see Appendix IB II for 
a more detailed discussion. Similarly to g^-type contributions, only a few of the N 2 LO terms 
do not cancel in the sum. Again, only those parts of the diagrams la, lb and mini-football 
that cannot be reduced to the topology of the football diagram in Fig. [31 give a non-zero 
contribution to the transition amplitude. The results are shown in Table [TTJ and the sum of 
these Qa contributions gives the following transition amplitude: 

*< L ° = % /<g^(M) {(2,+ - 2r_ - *0=£&oa • f 



We now turn to the contribution emerging from the diagrams of Fig. [3] with the c 2 and 
C3-vertices in the off-shell pion kinematics at nucleon 2. We obtain the following expression 
for the amplitude: 

x [dH c 3 (l-l) + (c 2 -g A /8m N )(v.l)(y.l) f 1 1 _ J 
7(27r) 4 (P-ml+iO)(P-ml + iO) I 2 2 J ' 

where the numbers in the curly bracket correspond to the individual contributions of the 
qa~ diagrams , as they appear in Fig. [3j in order. Again, while the individual diagrams do 
contribute at N 2 LO, their sum turns out to yield a vanishing result. We, therefore, conclude 
that there are no loop amplitudes oc q to the order we are working. 

Upon performing some simplifications, the total result for the ^-contribution to the 
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TABLE II: Interference pattern of s-wave contributions from the individual gA diagrams. The Table 
shows the contributions to the vertex e^T^Ag^ defined in Fig. [6] and Eq. (|16p . The contributions 
are given separately for different spin- momentum structures of the vertex A^, shown in the 
leftmost column. The notation for the isospin structures is defined in Eq. (|ll|) . 







Football 


Type la 




Type lb 


MiniFB 




Sum 


Order 


Si 


I 


-2r x 


-T+ + T_ + 


7"X 


T+-T- + T X 










NLO a , N 2 LO 


Si 


I 


-2r x 


-T+ + T- + 


T"X 


T+ - T- + T X 










NLO a , N 2 LO 


Si 


I v 


2r+ - 2r_ 


— T+ + T- + 


T"X 


-T + + T_ - T X 










N 2 LO 


Si 


l 


-2t + + 2t_ 


T+ — T— — 


Tx 


T+ - T_ + T X 










N 2 LO 


Si 


(Pl+Pl') 2^ 


2t x 


t+ - r_ - 


T"X 


-T+ + T_ - T X 










N 2 LO 


Si 


(Pl+Pl') 2^ 


2t x 


7^- — T— — 


Tx 


-T + + T- - T X 










N 2 LO 


Si 


-vl + iO 










2r+ - 2r_ + 2r x 





2r+ 


- 2t- + 2-r x 


N 2 LO 


Si 


f 2vq 
-vl + iO 





2t + - 2t_ - 


2r x 








2r+ 


- 2r_ - 2r x 


N 2 LO 


Si 


1 k'j-ml+iO 













8t x 




8r x 


N 2 LO 



"Notice that, in addition to the cancellation shown in the Table, in the case of NLO even the individual 
contributions to the corresponding spin-momentum structures turn out to vanish, see Appendix IB II for 
details. 



transition amplitude in Eq. ffTS]) can be brought into the form 



zM 



N 2 LO 



.9a[ ( w ,fdH v .(i + i) Sl .(l + l) 

7(2vr) 4 (P-ml + iO)(P-ml + iO)(-v ■ I + iO) 



-t x {v ■ q)(Si ■ h) 
+2r x (5 , 1 -A; 1 



dH 



v-(l + I) 



(2tt) 4 (p -ml + iO){P -ml + iO)(-v -1 + iO) 
I f dH v {l + l)q- (l + l) 



k\ -ml + iOJ (2vr) 4 (/2 - m 2 + l Q)(j2 - m l + iO) 



(20) 



The first term in Eq. fl20|) does not contribute at N 2 LO, since at this order the term f •(/+/) w 
2t> • / in the numerator cancels with the nucleon propagator — v • 1+iO. The resulting integral 
vanishes due to the symmetry of the integrand. Specifically, the integral is to be invariant 
under the shift of variables (/ — > —1,1 — > —I). Indeed, the denominator of this integrand 
is invariant under this transformation whereas the numerator changes its sign. Therefore, 
the first term in Eq. f j2"U|) is equal to zero. Finally, keeping only the lowest-order terms as 
appropriate at N 2 LO and using the expressions for the loop integrals outlined in Appendix IU| 
we arrive at the final result: 



%M. 



N 2 LO _ 9 A 



Hi 



— t x (v ■ q)(Si ■ hi) 



\j{kl) 



1 



18(4tt) 2 _ 



(21) 



where the UV-divergent integral J{k\) is defined in Eq. ffl~5]) . 



V. SUMMARY OF THE TWO-PION EXCHANGE DIAGRAMS 

Until now we have evaluated the expressions for the production operator assuming that 
the pion is produced from nucleon 1. We now add the contribution emerging from inter- 
changing the nucleon labels. We use the fact that in the center-of-mass system pi = —p 2 = p 
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and k\ = —k 2 + q and employ the approximate relation k\ ~ k\ with higher-order terms 
being ignored. Throughout, we also ignore operators leading to pion p-wave production. We 
then obtain from Eqs. (1211) and (1141) the following complete (i.e. symmetrized with respect 
to the nucleon labels) expressions: 



iM N 2 LO _ 9A (y ■ q) 



at 



f 5 



t x (S 1 + S 2 ) ■ h 



6 v lJ 18 



(47T) 



(22) 



9a 



Jit 



+r x (S 1 + S 2 )-k 1 



24 J(fclj + 9(47r)2 



(23) 



Employing dimensional regularization, d = 4 — e), the integral J(ki) entering the above 
expressions can be written in the form 



/r 



W) = - 



i 



i J (2tt)( 4 - £ ) [P -m 2 n + iO] [(I + h) 2 -ml + iO] 



-2L 



1 



(An) 2 

where the function FUx) is defined via 



(24) 



V4 - x - iO 
rAx) = 1= arctan 



x 



'x \ a/4 — x — iO 

and the UV divergency appears as a simple pole in the function L: 



(25) 



1 



(4tt) 



-- + - ( lE - 1 - log(47r)) 



(26) 



Note that both M„ N LO and M N 3 LO are proportional to the outgoing pion energy v ■ q ~ m^, 
i.e. both operator amplitudes vanish at threshold in the chiral limit. 



VI. REGULARIZATION PROCEDURE 

In MCS the loop diagrams which contribute to the renormalization of e.g. the nucleon 
mass m/v and the axial coupling constant gA do not involve large-momentum components. 
Consequently, these diagrams contribute in the MCS at order N 4 LO which is beyond the 
scope of the present work. For example, consider a LO rescattering diagram which in our 
naive counting is of order ^Jm^/m^. Including a pion loop in any of these diagram will 
require a renormalization any of the vertices in these LO diagrams, cf. e.g. the last three 
diagrams in Fig. HI This pion loop will increase the MCS order by factor (m^/m^) 2 as shown 
in, e.g. Ref. |2j, Table 11. At N 2 LO, we only have to consider the loop diagrams which are 
evaluated in this paper. The UV divergences appearing in the corresponding integrals are 
to be absorbed into LECs accompanying the 4iW amplitudes v4ct and £>ct introduced in 
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Sec. Ill Al The contributions of the loops to the amplitudes A and £>, see Eq. 
separated into singular and finite parts 



can be 



A 
B 



where 



(A 



singular 



+ A 



finite J 



singular 



+ B 



finite J 



^singular ^a(4"*) ~^ / ' ^singular 



9 i (^-1)(4.) 2 L. 



(27) 



(28) 



Here we have used that at threshold k\ = p and v ■ q = m v . Notice that the above de- 
composition into singular and finite pieces is, clearly, scheme dependent. Analogously, the 
amplitudes given by the 4Ntt Lagrangian contact terms, which are given in, e.g., Ref. ET 
are written as: 



A 



CT 



(47TA) a J? 



(A r c M + (4tt) 2 ^L), B C t 



m Tl 



(4vr/ 7r ) 2 / 7r 3 



(B r c M + (4tt) 2 /3 b L) . (29) 



The singular parts of the amplitudes in Eq. ( 12 9 p cancel the singularities of the amplitudes 
in Eq. (1271) . emerging from the loops. The resulting finite expressions for the scattering 
amplitudes are given in terms of the renormalized LECs of Ref. 



A 



m^ 



CT 



B 



CT 



(4tt/.) 2 £ 
m^ 



A r CT = -{d' x + 2ei - 2e 2 ; 
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^ T = -K + 2 ei ; 



(30) 



The magnitudes of the amplitudes *4. CT and £>£t can be estimated using the values of the 



LECs determined in Refs. [14J . I17L |37| where the short-ranged production mechanisms were 



assumed to originate from z-diagrams with a and u exchanges (see explicit expressions for 
these exchanges in Refs. 14, 29]). Given the estimates in Ref. [13], we find d[ + 2e\ — 



2e 2 ~ — 7.5/ f%m N and d[ + 2ei ~ — 3.5/ f%m N , and using m N ~ 47?/^, we obtain Aq T 
2m 7r /(m 2 v ^) and B r CT ~ lm^/{m 2 N f^), which results in A r CT ~ 2 and Bq t ~ 1. 

We take these numbers to set the scale for typical N 2 LO contributions. Therefore, these 
estimates allow us to infer the importance of the pion-nucleon loop contributions to the 
NN — > NNir reactions at threshold. In particular, we can compare this estimate with the 
finite parts of the loops given by Eqs. f[22]) and ([23]) (where v ■ p ~ m n <C |p|). 



"4finitc 
^finite 



9l 



9a 
' 6 



log 



ml 



2Fx 



-p- 



mi 



19 



g\ - 1 1 - log 



ml 



-2F X 



-V 



mi 



+ 3^~6 



(31) 



Choosing fi = Amf-^ with /„- = 92.4 MeV and qa = 1.32, we find -4. nn ite = —2.9 and 
-Sfinite = 1-4. We, therefore, conclude that contributions of the finite parts of the loops 
are comparable in size with A r CT and £> CT . This confirms our power counting and shows 
that pion loops contributions, not considered in previous analyses, are indeed significant. 
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One should, however, keep in mind that this result was obtained for the particular regu- 
larization scheme as explained above. In general, the finite parts of the loops ^4fi n ;te and 
Bfinite can be further decomposed into the short- and long-iaiage parts. The former one is 
just a (renormalization scheme dependent) constant to which all terms in Eq. (1311) but Fx 
contributes. On the other hand, the long-range part of the loops is scheme- independent. By 
expanding the function Fi(— p 2 /m 2 ), Eq. fl25|) . which is the only long-range piece in fl3Tj) . 
in the kinematical regime relevant for pion production, i.e. (p 2 /m 2 ) 3> 1, up to the terms 
at N 2 LO one obtains 

JSL = - 9 fio S (<) + o(< 



B£= = ^(^-illogl^l+Ol^). (32) 



V J \v 
9a /19 o A , / m 2 \ ( m 2 



12 V 4 



V J \v 



Numerical evaluation of these terms gives A^^^ = 2.2 and £>|°^f e = —1.5. The scheme- 
independent long range part of N 2 LO pion loops appears to be as large as the resulting 
short-range amplitudes, Aq T and Bq T , which are given by the meson-exchange mechanism, 
proposed in Refs. [(Jl! to resolve the discrepancy between phenomenological calculations 
and experimental data. Hence, the importance of the N 2 LO pion loop effects, not included 
in the previous studies, raises serious doubts on the physics interpretation behind the phe- 
nomenologically successful models of Refs. (6|-@] . 

In a subsequent work we will present results for N 2 LO loops including the Delta resonance 
as well as the convolution with proper nuclear wave functions. At that point a fit to the 
pion production data is possible and we can extract the strength of the counter terms from 
data. 



VII. SUMMARY AND DISCUSSION 



Chiral perturbation theory has been successfully applied in the past decades to describe 
low-energy dynamics of pions and nucleons. Application of this theoretical framework to pion 
production in nucleon-nucleon collisions is considerably more challenging due to the large 
three-momentum transfer involved in this reaction. The slower convergence of the chiral 
expansion for this reaction, i.e. the expansion in the parameter x ~ \Z m n/ m N defined in 
Eq. ((7]), provides a strong motivation for extending the calculations to higher orders. In 
this work we used the power counting scheme which properly accounts for the additional 
scale associated with the large momentum transfer, namely the momentum counting scheme 
(MCS), to classify various contributions to the NN — > NNn transition amplitudes according 
to their importance. We also evaluated all loop diagrams with pions and nucleons as the 
only explicit degrees of freedom up to and including N 2 LO. The considered loop diagrams 
can be divided into two groups according to the power of the nucleon axial-vector coupling 
constant: the ones linear with qa and the ones proportional to g A , see Fig. [3j We confirm 
the earlier findings that there are no NLO loop contributions to the threshold NN — > NNw 
reaction amplitudes. Our results, which are partially summarized in Tables |T] and [TT] and 
in Sees. |V] and IVI[ demonstrate that the MCS combined with the requirements of chiral 
symmetry (breaking) pattern of QCD lead to a high degree of cancellation among various 
N 2 LO contributions. In particular, all 1/m n- correct ions of the various diagrams cancel at 
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N 2 LO. We also show that the LECs q, i — 1 . . . 4, of do not contribute to the pion loops 
at this order. 

From Table [I] we see that only the cross-box diagram (diagram III) and the four-pion 
interaction diagram (diagram IV) contribute to the pion s-wave transition amplitude M^ 3 LO 

given in Eq. (I23"j) . The two cross-box diagrams contribute to both amplitudes, the isoscalar 
one A and the isovector one B, whereas diagram IV only contributes to B. Analogously, 
from Table HT1 one can deduce that the non- vanishing contributions to the amplitude M^ LO , 
Eq. fl22|) . originate from the double nN scattering diagrams of type la and lb and from the 
mini-football diagram. These diagrams however contribute only to the isovector amplitude 
B, as seen in Eq. fl22|) . Thus the only contribution from pion loops to the isoscalar amplitude, 
A, originates from the cross-box diagrams. 

The pattern of cancellations discussed above has important phenomenological implica- 
tions. In fact, none of the previous phenomenological investigations take into account either 
the cross box diagrams (type III) or the double scattering contributions (type I) which, as we 
find, contribute significantly to the production amplitude. In particular, the regularizat ion- 
scheme independent long-range contribution of the pion loops to A turns out to be com- 
parable in size with the short-range amplitudes emerging in phenomenological models of 
Refs. from heavy-meson z-diagrams which, in these studies, are advocated as the nec- 
essary mechanism to describe experimental data. Thus, our findings raise doubts on the 
role of the short-range physics in pion production as suggested in these phenomenological 
studies. We, however, refrain from making a more definite conclusions until the complete 
N 2 LO operator convoluted with the nucleon wave functions is confronted with experimental 
data 1 381. 



Meanwhile, within various meson-exchange approaches [ll|, [12|, |39|, the pion production 
is largely driven by tree-level pion rescattering off a nucleon with the irN — > irN amplitude 
being far off shell, see Fig. [71 The physics associated with irN scattering near threshold is 
normally parameterized in phenomenological calculations in terms of the a- and p-meson- 
exchange contributions. The scalar-isoscalar (a-type) irN interaction is relevant for the 
isoscalar production amplitude A while the isovector (p-type) tiN interaction contributes to 
the strength of B. The isoscalar nN scattering amplitude essentially vanishes on-shell, see 
Refs. [lCj for the most recent evaluation of the isoscalar tcN scattering length. Therefore, 



the mechanism of Refs. Il2l . [39| relies on the significance of the off-shell properties of 



the 7rN scattering amplitude. Our EFT consideration puts this mechanism into question. 
Pion rescattering via the phenomenological pion-nucleon transition amplitude can in chiral 
EFT be mapped onto pion rescattering (at tree level) via the low-energy constants q plus 
some contributions from pion loops. The tree-level piece oc q is, even in the off-shell (pion 
production) kinematics by far too small to explain the data for the neutral pion production 
[Til . [l5| . As far as the loop contributions are concerned, only diagrams IV and mini- football 
may be regarded as an analog of the corresponding phenomenological mechanism, as illus- 
trated in Fig. [7J However, after the cancellations, the only contribution that survives from 
diagram IV has an isovector structure as shown in Table [IJ Furthermore, the mini- football 
diagram gives an isovector contribution, see Table [Til Therefore, none of the pion loops can 
be mapped into the particular phenomenological mechanism in the isoscalar case. Thus, 
another lesson we learn from our work about the phenomenology of the neutral pion pro- 
duction near threshold is that the rescattering contribution with the isoscalar irN amplitude 
modeled phenomenologically by a a exchange should be very small. On the other hand, the 
rescattering mechanism with the isovector tiN amplitude, the Weinberg-Tomozawa term re- 
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FIG. 7: Diagram on the left-hand-side represents the phenomenological rescattering mechanism 
via the off shell 7r ^-amplitude. Diagrams on the right-hand-side are the only loop graphs at N 2 LO 
that can be interpreted as an analog of this phenomenological mechanism. For notation see Fig. [TJ 

lated to the p-meson exchange via Kawarabayashi-Suzuki-Riazuddin-Fayyazuddin relation 



4d . |41| , is potentially capable of resolving the discrepancy with the experimental data for 



charged pion production [2 
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Appendix A: The evaluation of the individual (/^-diagrams 

In this appendix we derive the NLO and N 2 LO expressions for individual two-pion ex- 
change diagrams shown in Fig. [3] under restriction that the outgoing pion is produced in 
s-wave. The kinematics is defined in Fig. 



1. Diagram II 



Diagram II shown in Fig. [3] is straightforward to evaluate. The operator A abc of Eq. (Q 
in this case arises from a three-pion one-nucleon vertex whose explicit form can be found in 
Eqs. (J5j), (J6j). The diagram II yields the contribution 

iM u = B 2 (l, I) r^r b |±J [rf 6 bc S l ■ (-1 + I) + r b 1 6 ac S 1 ■ (q + I) + r c l 6 ab S 1 ■ (q - I)} 

' ie abc [v ■ qS 1 ■ {-I -l)-v ISi -{l-q)+v- IS 1 ■ (q + I)} 



2m^ 

7^S 1 -(p 1 + p{) [r?5 bc v ■ (-Z + + r\ 5 ac v ■ (q + I) + r^5 ab v ■ (q - I)] 



19 



where Bi(l,l) is defined in Eq. ffTUj) . Contracting isospin indexes and ignoring all p-wave 
terms (oc S\ ■ q) and higher-order s-wave terms oc v ■ q/m^ ~ m^/rriN we find 

zM n = I—. I 5 2 (Z,ZW-(5i-Z)[4r + + 4r_-2r x ] + (5i-Z)[4r + + 4r_ + 2r x ] 



4/. 3 7 (27T) 

5 X -Z 



+ 



■u • Z + 1> • Z 



2rajv 2mjv 

5i-(pi+p0 



(2r+ - 2r_ 



2m, 



'/V 



[v ■ Z(4r+ + 4r_ - 2r x ) - v • Z(4r+ + 4r_ + 2r x )] 



where the integrand in the first line starts to contribute at NLO while that in the last two 
lines gives N 2 LO contribution. 



2. Diagram Ilia 

The crossed box diagram, Type Ilia, shown in Fig. |3] has a more complicated structure. 
In this diagram the operator consists of a ttN — > ttN scattering vertex, a nucleon 

propagator and a nNN-vertex. Again we only need to include the contributions from the 
leading and subleading chiral Lagrangian in the vertices. We also include the nucleon recoil 
correction in the nucleon propagator. This diagram gives the following expression 

d ^ d a a ~cjb 1 J>add (., n , „\ (Pi+Pi -l)-(l + Q) 



*M IIIa = J^L B 2 (l, I) rfr 2 6 ^- 2 s bad rt (v(l + q) 




9A ~ ci Sl -i- Sl ' {2pi ~ ' ,v - ). (Ai) 



We have ignored here the subleading Cj-contributions to the ttN — > ttN rescattering vertex 
since they are suppressed in the momentum counting scheme due to the negligible kinetic 
energy of the outgoing pion with q ~ (m^O). We will rewrite the ttN — > ttN vertex 
expression in the integrand above in a way similar to the rearrangement in Eq. (jHJ) 

„. (/+9 )_(a_+tf-o-( r +*>- 



2m, 



where we used that v ■ p[ = p[ ~ p{ 2 /2mjy. The first term on the second line of Eq. (|A2[) is 
identical to the nucleon propagator in Eq. f lAip and will give a factor of —1 when inserted 
into Eq. flAip . This factor of —1 together with the lowest-order contribution of the ttNN- 
vertex, Si ■ Z, give the NLO contribution of diagram Ilia. The last term in the second line of 
Eq. (IA20 contribute to an outgoing p-wave pion and is ignored in this paper. The 2q term in 
Eq. flA2p will contribute to the N 2 LO amplitude. We next use the relation Z = l + p\—p[—q 
in the 7rNN-vertex and in the nucleon propagator. We ignore p' w ~ go ~ m n contributions 
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and the recoil correction in the propagator which are of a higher order. Carrying out the 
isospin algebra we get: 

S ' • (■» + Pi) („ . f) _ f -^-O f 5l • i) | \2r + - r,] 



2rriN \—v-l + iO 

The first term in the curly bracket starts to contribute at NLO. The remaining three terms 
contribute to N 2 LO. 



3. Diagram Illb 

Diagram Illb (Fig. [3]) has a structure similar to diagram Ilia. We proceed along the same 
lines as for the two previous diagrams and obtain the following contribution: 

iM f d * 1 d n T\ ~cj> ( ,\9ah( q , Si • (2p[ - l)v ■ l \ 
%Mmh = /(2^ 52(M) T2T2i ~ 1] J l V 1 ^ ) 

i 1 ^ ( . ( _ r + ff) (pi + Pi - • H + g) 



Using the on-shell condition for the incoming nucleon with p 10 = pi 2 /2m/v, we rewrite the 
irN 7riV vertex in a way similar to what was done for diagram Ilia 

t,.(-r +g) -fo + fl- f >(- r+ fl = 

" /o " ~2^~ ) +2q ° 2^ • (A3) 

Using the relation 2q = 2v ■ q and keeping only terms appropriate at the order we are 
working we obtain: 

iM ™ = Wi /(^ B * f '{ (Sl -' ) + S (t "° 

- Sl t' +p;) (" • + ( w ■ '>}[-*-+ - *l 

2tyin \—v-l + iOJ J 

The first term in the curly bracket starts to contribute at NLO. The remaining three terms 
contribute to N 2 LO. 



4. Diagram IV 

Diagram IV (Fig. [3]) has an operator A a \ c containing a four-pion vertex, a pion propagator 

Q A. 

and one 7riViV-vertex. We keep the leading and next-to-leading order in the 7riViV-vertex 
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and obtain 



i Mi 



IV 



dH 



(27T) 



5 2 (/,Z) r 2 c r. 



c_6 / 9A 



fn 



Si ■ (pi +p[)v ■ ki 



2m, 



■N 



k[ — m 2 + iO 



x- 



f 2 

J IT 



[5 ab 5 cd [(I - qf - ml] + 5 ac 5 bd [(l + qf - m 2 ] + rt bc [(A* - g) 2 - m 2 ] } . 



The four-pion vertex is rewritten as a sum of six terms 



±{S*6« [(l-q) 2 -ml] + 

= A_{ 5 ab 5 cd [ l 2_ rn 2] +S ac 5 bd 
+S ab S cd ^_ 2l . g + g 2j + 5 ac 5 , 



(l + qf -ml 



72 2 

/ - m„ 



aczbd 



2l-q + q 7 



+ rt bc [(h - qf - ml] } 
+ 5 ad 6 bc [kf - ml] + 



(A4) 



The contributions from the first two terms on the r.h.s. of Eq. (IA4j) are of a higher order. The 
reason is that each term cancels a corresponding pion propagator in the operator ^(M)- 
When one pion propagator in I?2(M) is eliminated, the large momentum, like k\ or pi, of 
this reaction is no longer part of the loop integral which, consequently, only contributes at 
a higher order than what is considered in this paper. Keep in mind that v ■ k\, v ■ pi and 
v ■ P2 are all of the order m n , whereas, |fci| ~ p = ^/m^m^. The third term cancels the pion 
propagator k\ — m\ + iO and will contribute at NLO and higher order in our counting. The 
last three terms in Eq. (1A4j) start contributing from N 2 LO. 

Using k\ = I — I + q, dropping terms contributing to outgoing p-wave pions and carrying 
out the spin and isospin algebra we find: 



iMi 



fjA 



IV 



d 4 l 



4# J (2tt 

+ (5! ■ h) 



B 2 (l,l) 



(Si • - (Si • + 



Si ■ (gi + pj) 
2mjy 



-v ■ I + v ■ I 



6(r + + r_ 



k\ — m\ + iO 



(-8r x ) 



The first two terms in the first square bracket are NLO contributions. The remaining two 
terms are N 2 LO terms. 



5. Box diagram a 

In the expression for the Box a diagram (Fig. [3D we again rewrite the pion-nucleon rescat- 
tering vertex as a sum of two terms similar to what we did for the Type-Ill diagrams. One of 
the new terms will cancel nucleon propagator yielding an irreducible NLO contribution. In 
contrast to the derivation of the amplitude for the Type-Ill graphs, we here do not consider 
the contribution from the term with the (remaining) nucleon propagator since it is reducible 
and thus included in the initial iViV state interaction. Using again that the sum of the two 
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lowest orders contribute to the vertices, we obtain from the box a diagram: 



zM Box a = /77—T7 S 2 (Z, I) t 2 t 2 — £ r x \v ■ (-1 + q) 



(27r) 4 <v w , , 1 ^ v 

x f 5l . / - gi'(2pi + 0^^ (A5) 



2mjv 



To rewrite the expression in the pion-nucleon rescattering vertex we again use that that p[ 
is on-shell, i.e. p[ = p[ 2 /2rriN- The tcN — > ttN vertex is rewritten as 

v . { _ !+q) _(A±£ + ?)-(- f +d- 




- Mo + PM - ) + 2g ^ • 

The first term on the r.h.s. of the above expression is identical to the nucleon propagator 
and will give a factor of —1 when inserted into Eq. f lA5j) . The last term is a p-wave pion 
contribution and is ignored. Also the 2g -term does not need to be taken into account as 
it corresponds to a reducible contribution. Using 2pi + I = I + (j>i + p{) + q, ignoring the 
p-wave pion terms, and evaluating the spin and isospin structures, we find: 

rirred. 9A f d 4 l ~ f S l'h.. A , S l ' fa + Pi) 



6. Box diagram b 

The Box b diagram given in Fig. [3] is very similar to the Box a diagram and the evaluation 
procedure is similar. We consider again only the irreducible contribution. The diagram gives: 



x ? (v ■(( + ,)- (PL±S±M±«1 I . (A6 ) 



2mjv 

Again, rewriting the pion-nucleon rescattering vertex using that p± is on shell, pio = pi/2rriN 
leads to : 

„. (J+9) _®±s±6j±a = 

+ /o " J + 2go 2^ • (A7) 

The first factor on the r.h.s. of the Eq. (IA7I) . when coupled with the nucleon propagator in 
Eq. (IA6I) . yields a factor of 1 while the 2g -term in Eq. (|A7[) produces a reducible contribution 
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included in the final NN state interaction. Using the relation 2p[ + I = I + (pi + p[) — q, 
ignoring terms leading to outgoing p-wave pions and carrying out the spin and isospin algebra 
leads to: 

Like the final expression for the Box a diagram, the first term starts at NLO and the next 
two terms are the N 2 LO contributions to the amplitude. 



Appendix B: The evaluation of the individual ^-diagrams 

In this appendix we derive the expression for two-pion exchange diagram linear in for 
s-wave pions produced. The final expressions for the diagrams contain N 2 LO contributions. 
The kinematics is defined in Fig. [6j 



1. The Football diagram 



The two pion propagators are tied together in pion-nucleon scattering vertices at both 
nucleons. Since this loop diagram involve just pion propagators, we have an extra symmetry 
factor 1/2 associated with the boson loop. The football diagram shown in Fig. [3] gives the 
following expression: 



iM F 



x | [r^6 bc S 1 ■(-/ + /) + r b 1 6 ac S 1 ■ (q + I) + r c l 5 ah S 1 ■ (q - I)} 
' ie abc [v ■ qSi ■ (-1 -l)-v- lS t -(l-q)+v- IS 1 ■ (q + I)} 



2m N 
1 

2mjy 



Si ■ (Pi + piM5 bc v ■ (-1 + + T b 6 ac v -(q + l) + r{b ab v • (g - 0U- 



After performing some spin and isospin algebra, dropping terms corresponding to the out- 
going p-wave pion and/or higher-order corrections we obtain 



(Bl) 



+ 



2rriN 2m,N 



(2r + - 2r_) 



Si ■ {pi+pD 



2m, 



TV 



V ■ I + V ■ I 



(2r x 



The obtained result requires some clarification. Looking naively at the first line in Eq. flBT 
one may conjecture that this diagram starts to contribute already at NLO. Indeed, assuminj 
l ~ \l\ ~ p, the dimensional analysis gives 



P 



p 



1 p 2 



f 3 f 2 p 3 (47r) 2 £ 2 



N 
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where the three terms on the l.h.s. stand for the 37riViV-vertex, the estimate of D(l,l) 
as follows from Eq. (|T7|) and the integral measure, in order. Above we also used that 
(4irf n ) 2 ~ m^. On the other hand, a more careful analysis shows that the first line of the 
integral in Eq. (IBip . which appears at NLO, is 



• w 9a f dH 2l 

iMy = r x 



16 ft J (2tt)* (] 2 



(l 2 -l 2 + i0)(l 2 -l 2 + i0) 



(S 1 ■ I) + (Si • I) 



where we have used that Iq ~ \l\ ~ p 3> m^ to drop all sub leading contributions including 
the 1/mjv terms in the curly bracket of the integrand. This last integral, however, vanishes 
when integrating over Iq because the numerator of the integrand is an odd function of Iq 
whereas the denominator is an even one. The next-higher order contributions in Eq. (IBip 
do not vanish. They scale as m n /p and p/m^ compared to NLO and thus emerge at N 2 LO. 
Following the same lines, one can show that also the other diagrams of p^-topology start to 
contribute at N 2 LO. 



2. Diagram la 

The double-scattering diagram la shown in Fig. [3] gives the following expression 
dH n n h^cbyv ( 1 \ badd L /, , x (Pi +P{ -l)-{l + Q) 



lMl * = D ^ ~ l)eChVrl (4) £badr " [ V -V + ^- 



2m, 



N 



pio _ l o - ~2^T + M \ 

In the ttN — > ttN rescattering vertex (off nucleon 1) we included the leading WT vertex 
contribution together with its recoil correction. However, we dropped the subleading q- 
terms in this vertex since they are of higher order (see the discussion in the end of Sec. IIIBl) . 
The tcN — > tcN vertex expression is rewritten the same way as for diagram Ilia, shown in 
Eq. flX2j). 

Using that p[ = p[ 2 /2m N , collecting the spin structures and performing the isospin 
algebra we get: 



Si ■ (pi+pi) 



2m, 



v ■ I + v ■ I 



Av ■ q 



-v -1 + iO 



(sW~)}(t + -t_-t x ). 



This amplitude starts to contribute at N 2 LO, see discussion in Appendix IB II for more 
details. 
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3. Diagram lb 

The double-scattering diagram lb in Fig. [3] gives an initial expression: 
iM lh = 



dH 



(2*) 



d 2 (m> c6 M(-i)^i [Si-i- 



J5r 



Si • (2p{ - Z)v • Z 



r 



(pi+Pi-Z)-(-Z + g) 



2m 



■JV 



Again, the 7rA^ — > ttN vertex is rewritten as sum of two terms as for diagram Illb, see 
Eq. ( ]A3j) . We follow the simplifications discussed for diagram Illb, and use that p w = 

Using again Z — » — Z etc. to simplify the integrals containing the function D2(l, I), collecting 
spin structures and performing the isospin algebra we get: 



fl4 



lb 



dH 



Iff J (2vr) 4 
Si ■ (pi+p{) 



D 2 {l,l)\(S 1 -l) + (S 1 -l) 



2?% 



u-Z + 



5i-Z 


2mjv 



2m, 



t> ■ Z + 1> • Z 



+ 



4-y • g 



-•u • Z + z'O 



(SW) [ r+ -r_ + r x 



This amplitude also starts to contribute at N 2 LO, see comment in Appendix IB II for more 
details. 



4. The Mini-Football diagram 



The contribution of the mini-football diagram in Fig. [3] can be written as 

1 f dH „ r, „»,„. „ / Qa\ j(„ , Si ■ (pi + p[)v ■ h\ / 



mF 



2 7(2tt) / 



D 2 (l,l) s^ri [ g 4 K (S 1 -k 1 



2m N 



k\ — m\ + iO 



2 2 



x — {<r fe <5 cd [(Z - g) 2 - ml] + <5 ac 5 M [(Z + q) 
The four-pion vertex can be rewritten as a sum of six terms 



+ [(fc - g) 2 - m 2 ] ) . 



■1 [(Z - qf - ml] + ^ [(Z + g) 2 - m 2 



+S ab S cd f_ 2 / . q + q *] + gacg. 



I — mi 



aczbd 



2l-q + q 2 



+ 5 ad 5 bc [(h - g) 2 - ml] } 
+ 5 ad 5 bc [k\ - ml] 



Following the arguments outlined in the derivation of the contribution from diagram IV, see 
the discussion below Eq. (1A4|) . most terms either contribute to outgoing p-wave pions or 
higher orders in the chiral expansion. The final result reads: 



This amplitude starts to contribute at N 2 LO. 



g-(Z + Z~) 
k\ — m 2 + iO 



:8r v ). 
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Appendix C: Expressions for loop-integrals 



In this appendix we provide expressions for loop integrals required to calculate the tran- 
sition amplitude at N 2 LO. Using dimensional regularization and integration procedure de- 
scribed in Appendix E of Ref . [42] , we obtained the following results: 



dH 



v{l + l)S 1 ■ (I + I) 



i J (2tt) 4 (p -ml + iO)(P -ml + iO)(-v ■ I + iO) 
1 f dH v-(l + I) 



% J (2tt) 4 (/2 _ m 2 + iQ )(/2 _ m 2 + iQ )(_ v .l + i 0) 

dH v .(i + i) q .(i + i) 



■2J(kl 



i J (2vr) 4 (P -ml + iO){P - ml + iO) 



2kf(vq) 



6 J(A:?) 9 (Any 



(CI) 

(C2) 
(C3) 



dH 



(S 2 ■ I) (S 2 ■l)S 1 -(l + I) 



(2tt) 4 (p - ml + iO)(P -ml + iO)(-v ■ I + iO) 2 



i r dH 



(S 2 ■ I) (S 2 ■ I) 



i J (2tt) 4 (P -ml + iO)(P -ml + iO)(-v ■ I + tOf 



-te^k lfl S ll/Va S 2 pJ(k 2 ), (C4) 



\m) - ^ cos) 



dH 



(S 2 -l)(S 2 -l)q-(l + l) 



% J (2tt) 4 (p -ml + iO)(P -ml + iO)(-v -1 + iO) 



v ■ q 



12 J(A; ' )+ 18(4tt)2 



(C6) 



where integral J(k\) is given by Eq. f lTo^) . and only the leading loop contributions for the 
s-wave pion in MCS are kept. 
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